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ABSTRACT 


Moisture  data  from  radar  and  rawinsonde  observations  during  three  lake  effect  snow  events  are 
analyzed  to  determine  entrainment  rates.  Gathered  during  the  winter  of  1990,  the  data  is  a  sub¬ 
set  from  the  Lake  Ontario  Winter  Storms  (LOWS)  Project,  which  deployed  a  mesoscale  network 
of  sensors.  Doppler  wind  profiler  signal-to-noise  ratio  (SNR)  data  is  used  to  derive  humidity 
structure  function  parameter  (C?2)  time-height  series  analysis,  which  is  then  compared  to 
rawinsonde  specific  humidity  (q)  plots.  Visual  comparison  of  log(C?2)  and  q  analysis  indicated  a 
strongly  positive  correlation.  Radar-derived  humidity  analysis  is  used  to  estimate  the  depth  of  the 
Type  I,  cloud-topped  boundary  layer  (CTBL),  which  corresponded  well  with  results  from  LOWS 
rawinsonde  data.  Calculations  of  the  contribution  of  (C?2)  to  the  refractive  index  structure 
parameter  (C„2)  showed  the  humidity  correction  factor  (a2)  to  range  from  1.03  to  1.04  within  the 
CTBL,  consistent  with  previous  findings  for  Type  II  CTBLs.  A  comparison  of  entrainment 
rates,  computed  via  two  different  methods,  were  in  agreement. 
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1 .  Introduction 


Much  of  the  theory  behind  radar  observations  in  clear  air  comes  from  turbulence 
studies.  Researchers  in  the  1940s  through  the  early  1960s  laid  the  foundation  for  the  stoichaotic 
description  of  the  atmospheric  turbulence.  From  the  result,  clear  air  radar  observation  theory 
was  refined  in  the  late  1960s  through  the  early  1990s.  Operational  applications  rapidly  followed 
theoretical  developments.  It  was  discovered  early  that  the  turbulent  mixing  of  the  atmosphere 
causes  gradients  of  temperature,  moisture,  and  pressure,  and  hence  gradient  of  the  refractive 
index.  Radars  of  sufficient  wavelength  are  capable  of  sensing  these  refractive  index 
inhomogeneities,  as  they  are  almost  always  present  within  the  turbulent  atmosphere  below  the 
stratopause. 

This  thesis  presents  a  study  of  the  moisture  profiles  within  a  Type  I  cloud-topped 
boundary  layer  (CTBL)  using  wind  profiler  and  rawinsonde  data.  It  closely  mirrors  the  study  of 
the  moisture  structure  of  a  Type  II  CTBL  conducted  by  White  et  al.  (1991)  off  Southern 
California.  A  Type  II  CTBL  is  defined  as  a  convective  boundary  layer  primarily  driven  by 
radiative  cooling  at  cloud  top  (with  nearly  thermal  equilibrium  at  the  sea  surface).  The  Type  I 
CTBL  differs  in  that  it  is  driven  largely  by  heating  from  below  as  cold  air  moves  over  relatively 
warm  water  (Agee,  1990).  The  data  utilized  comes  from  the  Lake  Ontario  Winter  Storms 
(LOWS)  Project,  a  detailed  investigation  of  mesoscale  motions  within  Lake  Ontario  lake-effect 
snowstorms  (LESS)  during  the  winter  of  1990  (Reinking  et  al.,  1993).  In  the  LESS 
environment,  surface  fluxes  of  moisture  and  temperature  drive  convection  that  lead  to  localized 
snowfall  events.  The  LOWS  data  set  includes  in-situ  measurements  of  the  planetary  boundary 
layer  using  rawinsondes,  instrument  towers,  and  surface  observations  to  increase  knowledge  of 
dynamic  motions  and  meteorological  variables  within  the  mesoscale  environment.  Doppler  radar 
wind  proflilers  were  also  key  remote-sensing  tools,  having  the  advantage  of  increased  spatial  and 


1 


temporal  sensing.  From  the  radar  reflectivity  and  spectrum  width  data,  important  meteorological 
information  can  be  derived.  Reflectivity  data  caused  by  refractive  index  irregularities  in  clear  air 
(Bragg  scattering)  will  be  used  in  this  study  to  derive  entrainment  rates.  These  rates  will  be 
analyzed  in  regards  to  boundary  layer  growth  at  one  wind  profiler  site.  Comparisons  will  be 
made  with  the  results  of  Type  II CTBL  studies  and  with  previous  conclusions  from  the  1990 
LOWS  study. 

2.  Data 

Lake  Ontario  LESS  events  are  unique.  There  is  a  well  defined  snowbelt  as  a  result  of 
the  lake's  length  and  orientation  to  the  prevailing  winds,  and  due  to  the  orographic  enhancement 
effects  of  rising  terrain  on  the  leeward  (in  this  case  eastern)  portion  of  the  lake.  Single  band 
snow  storms  frequent  Lake  Ontario  and  Lake  Erie,  remaining  infrequent  events  on  the  other 
Great  Lakes.  However,  Lake  Ontario  is  much  deeper  than  Lake  Erie,  helping  to  maintain  large 
areas  of  ice-free  lake  surface  and  potential  for  heat  and  moisture  fluxes  into  the  boundary  layer 
throughout  most  winters.  On  the  other  hand,  Lake  Erie's  LESS  season  is  considerably  shorter,  as 
ice  often  covers  most  of  the  lake  surface  by  mid-winter. 

2.1  Lake-effect  Snowstorms 

Primarily  a  mesoscale  phenomena,  lake-effect  storms  are  however  a  complex 
interaction  of  synoptic  scale,  mesoscale,  and  microscale.  Lake-effect  snow  events  develop  when 
the  atmosphere  is  destabilized  as  post-frontal  Arctic  or  polar  air  passes  over  a  relatively  warm 
lake  surface.  This  causes  a  positive  flux  of  heat  and  moisture  into  the  atmosphere  producing 
convection  and  thickening  of  the  boundary  layer  through  entrainment.  Single  or  multiple  parallel 
bands  of  precipitation  form  which  can  produce  localized  heavy  snowfalls.  Wind  direction  plays 
an  important  role  in  determining  band  types  and  location.  Advection  length-wise  along  the  lake 
tends  to  produce  intense  single  bands,  while  all  other  wind  directions  produce  weaker  multiple 
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bands.  Additionally,  orography  serves  to  enhance  vertical  velocities  and  increase  snowfall 
amounts.  The  primary  and  secondary  factors  which  influence  lake-effect  snow  development  are 
listed  in  Table  1. 


;  Secondary  Factors  Related  to  Primary 

Strong  flow  of  Arctic  air 
across  a  relatively  warm  lake 

850  mb  temperature  of  -10°  i  -  850  mb  temperature  of  -12° 
for  lake  water  ^  2°  1  for  lake  water  <  2° 

Relatively  long  fetch  across 
lake 

40  miles  with  PVA  and  850  .  - 100  miles  otherwise 
temperature  £  -5°  . 

Positive  vorticity  advection 

NVA  and  rising  surface  , 

pressure  are  deterrents  to  ! 

snow  production  , 

Frictional  and  thermal 
convergence 

1 

1 

i 

1 

i 

||  Latent  heat  release 

1 

i 

Table  1.  Primary  LESS  and  secondary  related  factors.  Adapted  from  Penc  (1995). 

2.1.1  Lake-effect  Storm  Types 


Single  or  double  banded  storms  form  parallel  to  the  winds  aloft  and  aligned  along  the 
long  axis  of  the  lake  due  to  organized  convergent  flows  which  serve  to  force  convection  into  the 
center  of  the  lake.  Typically,  storm  dimensions  are  2-20  km  wide  and  50-200  km  long,  with 
cloud  tops  to  4  km  (Peace  and  Sykes,  1996).  The  single  banded  storm  is  the  most  intense  of  lake 
effect  snow  storms,  capable  of  producing  heavy  localized  snowfalls. 

Multiple  banded  storms  occur  when  the  winds  are  aligned  along  the  short  axis  of  the 
lake.  Normally  these  storms  are  5-80  km  long  with  aspect  ratios  ranging  from  5:1  to  12:1  (Peace 
and  Sykes,  1996).  The  boundary  layer  remains  relative  shallow  with  1-2  km  cloud  tops.  The 
bands  drift  with  the  wind  at  small  angles  to  the  mean  wind,  generating  widespread,  but  light, 
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precipitation. 

2.1.2  Thermal  Instability  Indicators 


Operationally,  the  Collier  index  is  used  to  determine  potential  for  and  severity  of  lake 
effect  snow  events,  measuring  the  degree  of  thermal  instability  between  the  lake  surface  and  the 
lower  troposphere  (Niziol,  1987).  Additionally,  it  has  been  noted  that  a  deepening  of  the 
boundary  layer  is  almost  always  associated  with  the  events  (Penc,  1995).  Boundary  layer  depths 
of  3-4  km  are  seen  with  the  most  severe  storms  (Niziol,  1987).  Numerical  modeling  studies 
have  shown  heating  and  moistening  of  the  atmosphere  are  responsible  for  increasing  boundary 
layer  depth  by  hundreds  of  meters  (Lavoie,  1972).  During  LOWS,  changes  in  the  height,  as 
well  as  the  strength  of  the  capping  inversion  were  indicated  to  be  the  most  important  factors  in 
lake  effect  snow  band  intensity.  LOWS  researchers  found  the  intensity  of  the  convection  was 
more  closely  related  to  the  depth  of  the  convective  mixed  layer  than  to  the  degree  of  lower 
atmospheric  thermal  instability.  The  greatest  snowfalls  occurred  during  single  band  storms  on 
12  and  14  January  1990,  where  CBL  heights  ranged  between  2.5  and  4  km  (Penc,  1995). 

2.2  Lake  Ontario  Winter  Storms  (LOWS') 

The  Lake  Ontario  Winter  Storms  project  was  conducted  from  5  January  to  1  March 
1990  to  demonstrate  and  evaluate  mesoscale  observation  and  forecasting  techniques  applicable 
to  location-specific  lake-effect  storms  and  freezing  rain.  In  order  to  resolve  mesoscale  features, 
an  array  of  six  specialized  remote  sensors  (a  dual-polarization  short  wavelength  radar, 
microwave  radiometer,  radio  acoustic  sounding  system,  and  three  wind  profilers)  were  deployed, 
as  well  as  a  full  compliment  of  project-specific  and  standard  meteorological  observing  systems. 
Researcher  from  various  universities,  government  agencies,  and  private  consulting  firms 
participated  in  the  project  as  summarized  in  Reinking  et.  al.  (1993). 
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2.2.1  Meteorological  Conditions 


Due  to  wanner  than  normal  air  temperatures  during  LOWS,  occurrences  of  lake-effect 
snow  events  were  below  normal.  To  establish  the  scope  of  this  study,  lake-effect  snow  data  sets 
were  analyzed  for  the  single  band  storms  of  11-13  January  and  14  January,  and  the  multiple 
band  storm  events  of  25-26  February  and  27  February  - 1  March.  This  data  encompassed  the 
four  "marginal"  or  greater  episodes,  as  described  by  Penc  (1995),  and  included  sub-hourly  radar 
and  rawinsonde  data  collected  during  LOWS  intensive  observation  periods  (IOPs). 

2.2.2  11-14  January  1990 

After  the  passage  of  an  occluded  front  between  2100  and  2200  UTC  on  1 1  January,  the 
synoptic  situation  was  favorable  for  development  of  lake-effect  snow.  During  12  January,  an 
intense,  980  mb  surface  low,  approximately  500  km  north  of  Lake  Ontario,  tracked  east 
producing  west  to  southwest  flow  along  the  major  axis  of  the  lake  from  0000  to  1 800  UTC. 

Aloft  at  500  mb,  a  broad  trough  also  produced  westerly  flow.  Surface  air  temperatures  ranged 
from  -1°C  to  -3°  with  the  Collier  Index  remaining  in  the  conditional  category.  Deep  cold  air 
advection  along  with  moderated  positive  vorticity  advection  (PVA)  produced  destabilization  that 
supported  development  of  a  stationary  single  snowband  over  the  lake.  Four  snowbands  formed 
over  the  center  portion  of  the  lake  beginning  around  0130  UTC,  with  the  final  snowband 
dissipating  around  1700  UTC.  The  snowbands  formed  over  the  warmest  part  of  the  lake  (near 
4°)  and  moved  inland  over  the  eastern  shores  of  the  lake  producing  up  to  80  cm  snowfall 
amounts. 

At  0000  UTC  on  14  January,  the  surface  low  had  moved  to  a  position  off  the  Canadian 
maritimes  and  was  deepening,  while  a  significant  shortwave  trough  stretched  from  James  Bay 
toward  Minnesota.  Surface  air  temperatures  ranged  from  -13°  northwest  of  the  lake  to  -6°  to  -9° 
southeast  of  the  lake.  The  Collier  index  remained  in  the  conditional  category,  but  approached 
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moderate.  The  approaching  shortwave  trough  and  broadening  longwave  trough  led  to  increased 
westerly  flow  and  longer  fetch  over  Lake  Ontario.  Snowfall  with  this  event  was  confined  to  the 
southeast  shore  of  the  lake  with  accumulations  up  to  15  cm.  Visibilities  dropped  to  zero  at 
Oswego,  NY  at  0600  UTC,  with  occasional  thunder  and  lightning. 

Of  note  is  the  possibility  that  some  of  the  snowfall  during  these  two  events  may  be 
attributed  to  synoptic  forcings,  as  well  as  to  lake-effect  snow  conditions.  The  resulting  snowfall 
may  be  referred  to  as  lake-enhanced  snow.  However,  for  the  scope  of  this  study,  only  the  lake- 
effect  factors  will  be  considered. 

2.2.3  25  -26  February  1990 

Multiple  lake-effect  snowbands  occurred  after  a  west-east  oriented  surface  trough 
passed  over  the  southern  portion  of  the  lake  approximately  1200  UTC  on  25  February.  Snow 
fell  over  the  entire  region  just  prior  to  the  trough  passage,  which  was  then  followed  by  light  lake- 
effect  snowfall  for  a  12  hour  period  over  the  southeastern  portion  of  the  lake.  A  tight  pressure 
gradient  providing  strong  northwest  winds  after  the  passage  of  the  trough  advected  in  very  cold 
surface  temperatures  (-20  to  -26°C).  Despite  a  short  fetch,  thermal  instability  was  great  enough 
produce  moderate  to  extreme  Collier  indices.  By  1500  UTC,  multiple  north-south  bands  were 
seen  over  the  southeastern  portion  of  the  lake  and  with  a  moderate  burst  of  snowfall  at  the  North 
Rose  radar  site.  As  the  boundary  layer  winds  continued  backing,  snow  band  orientation  changed 
through  northwest-southeast,  and  eventually  northeast-southwest  by  0300  UTC  on  26  February. 
Strong  surface  to  850  mb  ridging  was  evident  by  0000  UTC  on  26  February  causing  rapid 
slackening  in  the  pressure  gradient  and  decreasing  cold  air  advection.  By  1200  UTC,  only  a 
solitary  weak  dissipating  snow  band  remained.  In  general,  the  lake-effect  snowfall  totals  were 
very  light  in  this  short-lived  event.  Only  4  cm  accumulated  at  North  Rose,  directly  attributed  to 
lake-effects. 
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2.2.4  27  February  - 1  March  1990 


A  cold  front  approached  on  27  February  and  passed  over  the  lake  around  0000  UTC  on 
28  February,  after  which  multiple  snowbands  developed.  At  1200  UTC  on  28  February  there 
was  high  pressure  centered  over  the  upper  Midwest  producing  northwesterly  flow  over  the  Great 
Lakes  region.  Surface  temperatures  ranged  from  -12°  north  of  the  lake  to  -8°  southeast  of  the 
lake  while  850  temperatures  were  -17°  to  -18°,  leading  to  a  moderate  Collier  index.  With  a 
500mb  trough  providing  westerly  flow  aloft,  little  directional  shear  was  present  in  the 
troposphere.  By  1400  UTC  on  28  February,  radar  showed  multiple  snowbands  extending  inland 
on  the  southeastern  shore  of  Lake  Ontario.  By  1700  UTC  radar  echoes  ceased.  More  weak 
banded  echoes  resumed  at  2200  UTC,  then  dissipated  by  0400  UTC  on  01  March.  Snowfall 
along  the  southeastern  portion  of  the  lake  was  >  5cm  with  a  maximum  of  15cm  (Penc,  1995). 

2.3  Radar  Specifications  and  Siting 

The  radar  data  analyzed  was  retrieved  at  North  Rose,  NY,  on  the  southeast  comer  of 
Lake  Ontario  approximately  5  km  inland  from  the  lake  shore.  From  this  location,  measurements 
could  be  taken  within  and  near  single  and  multiple  banded  snow  events.  Multiple  band  snow 
events  frequent  the  North  Rose  site,  whereas  single  band  events  often  occur  farther  north, 
allowing  the  North  Rose  site  to  provide  measurements  of  parameters  associated  land-to-lake 
convergence  areas.  The  data  was  recorded  on  magnetic  media  by  the  operators  of  the  404  MHz 
UHF  Doppler  wind  profiler  provided  by  Penn  State  University.  This  profiler  used  a  three  beam 
system  with  two  orthogonal  off  zenith  and  one  vertical.  Data  archived  included  returned  signal 
power,  wind  velocity,  and  wind  spectral  width,  which  are  the  zeroth,  first,  and  second  spectral 
moments,  respectively.  Measurements  were  made  in  high  mode  (east,  north  beams)  followed  by 
low  mode  (vertical,  east,  and  north  beam)  every  3  1/2  minutes.  During  intensive  observation 
periods,  the  sub-hourly  information  was  recorded  for  all  range  gates  (Penc,  1992).  For  the  scope 
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of  this  study,  the  sub-hourly  signal  to  noise  ratio  (SNR)  data  from  the  low  mode  vertical  beam 


was  used. 

2.3.1  Data  Processing  and  Quality  Control 


In  order  to  accurately  calculate  atmospheric  variables  from  radar  reflectivity  data,  the 


Radar  Characteristics 

Frequency 

404.37  MHz 

Wavelength 

0.742m 

Maximum  bandwidth 

2  MHz 

Peak  power 

9.0  kW 

Duty  cycle 

<25% 

Antenna  Aperture 

47.0  m2 

Antenna  pointing 

zenith,  14.5°off-zenith  to  north  and  east 

Antenna  type 

Fixed-phased  array  of  collinear,  coaxial 

dipoles 

Antenna  element  spacing 

0.3125  m 

One-way  beamwidth 

-7° 

System  noise  temperature 

~ 100°K 

Operating  Parameters 

Vertical  Data  Processing 

Pulse  width  (t) 

1  |TS 

Pulse  repetition  period 

100  ps 

Average  power 

80  W 

Time  domain  averaging  (nc)  200  pulses 

Spectral  averaging  (nj 

8  spectra 

Maximum  radial  velocity 

±  9.27  ms-1 

Spectral  resolution  (IC3) 

0.290  ms'1 

Height  Sampling 

First  height 

0.40  km 

Height  spacing 

100  m 

Number  of  heights 

22 

Table  2.  Radar  characteristics  and  operating  parameters  for  the  Penn  State  404  MHz  wind  profiler 
deployed  at  North  Rose,  NY.  (Penc,  1992) 

radar  must  be  properly  calibrated.  The  Penn  State  404  MHz  wind  profiler  was  calibrated  at  Otis 
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Air  Force  Base  in  southeastern  Massachusettes,  1 1  months  before  the  LOWS  project  began. 
Calibration  used  data  from  an  aircraft-mounted  Lyman- a  hygrometric  sensor  that  was  compared 
with  radar-deduced  humidity  profiles.  This  radar  is  the  same  as  used  by  White  et  al  (1991)  in  the 
study  of  humidity  profiles  in  a  Type  II CTBL. 

Signal  processing  is  largely  based  on  techniques  developed  at  the  NOAA  Wave 
Propagation  Laboratory  (Staunch  et  al,  1984).  Before  archival,  the  raw  data  is  filtered  using  a 
series  of  steps.  First  is  a  coherent  integration  or  time  domain  averaging,  which  reduces  the 
amplitude  of  components  outside  the  maximum  frequency  measured.  Next  is  a  spectral  or 
incoherent  averaging  which  improves  the  detection  of  the  spectral  peak  by  smoothing  out  the 
noise  floor.  Then,  a  low-pass  filter  is  employed  which  rejects  data  below  -18  dB.  In  addition  to 
the  aforementioned  filtering,  data  quality  control  was  achieved  by  periodic  monitoring  of  the 
SNR  and  comparison  with  collocated  rawinsondes. 

3.  Remote  Sensing  Theory 

Radar  returns  in  clear  air  result  from  inhomogeneities  in  the  refractive  index  caused  by 
turbulence.  The  intensity  of  the  return  radar  caused  by  the  inhomogeneities  can  be  used  to  derive 
temperature  and  humidity  values.  In  the  study  of  propagation  of  electromagnetic  energy  through 
a  turbulent  atmosphere,  the  refractive  index  structure  parameter  ( Cn2)  is  the  key  to  analyzing 
returned  radar  backscatter  power  (Pr)  in  terms  of  the  refractive  index  (n),  temperature  (7),  and 
moisture  (specific  humidity,  q )  analysis.  It  has  been  shown  that  the  returned  backscatter  power 
received  by  the  radar  is  proportional  to  the  refractive-index  structure  parameter  (Van  Zandt  et  al, 
1978).  The  theory  behind  this  remote  sensing  is  presented. 

3.1  Turbulence  parameters 

Following  Fairall  (1991),  consider  some  atmospheric  variable  X  (r)  as  a  function  of 
location  r.  The  structure  function  DJd)  of  X  is  defined  by 
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DJd)  =  [X  (r  +d)  -  X  (r)]2  (3.1) 

where  d  denotes  a  separation  distance  and  the  overbar  denotes  an  ensemble  average.  Assuming 
isotropic  turbulent  fluctuations  within  the  inertial  subrange,  the  structure  function  parameter  can 
also  be  expressed  in  the  2/3-law  form  as  (Tatarski,  1961) 

C2  =  Dx(d)  /  (F3  (3.2) 

For  the  refractive  index  (n) 

Dn(d)  =  [n'(r  +d)  -  n'  (r)]2  =  C2  dm  (3.3) 
where  n'  indicates  perturbations  in  the  spatial  mean  of  the  refractive  index. 

3.2  The  Radar  Equation 

The  radar  equation  can  be  expressed  as  a  relation  of  backscattered  power  at  the  radar  antenna  Pr 
to  the  average  reflectivity  per  unit  volume  t]  of  the  target 

Pr=  P,  G2A2ff  ht]/(  1024 7^ (ln2)  r2)  (3.4) 

where  P,  is  the  radar  transmitted  power,  Gis  antenna  gain,  A  is  the  radar  wavelength,  0  is  the 
radar  beamwidth  (full  width  at  half  maximum)  in  radians,  h  is  the  transmitted  pulse  length,  r  is 
the  range  to  the  target.  From  the  radar  equation  a  value  of  t]  can  be  obtained,  since  Pr  can  be 
measured  and  the  remainder  of  the  radar  equation  parameters  known.  Again,  assuming  the 
refractive  index  fluctuations  are  a  result  of  isotropic  turbulence  in  the  inertial  subrange  and  of 
length  scale  within  At 2,  the  reflectivity  of  clear  air  can  be  expressed  as  (Ottersten,  1969) 

/?  =  0.38C2A'1/3  (3.5) 

From  this  expression,  it  can  be  observed  that  for  radar  wavelengths  of  10  cm  or  more,  the 
structure  constant  C2  is  directly  proportional  to  the  radar  reflectivity  r).  Radars  of  shorter 
wavelengths  are  sensitive  to  smaller  scales,  and  therefore  may  be  observing  contributions  of 
dissipation  by  molecular  diffusion,  depending  on  the  energy  dissipation  rate  (Ottersten,  1969). 


10 


3.3  Humidity  Structure  Function  Parameter  (r  2\ 

The  radio  refractive  index  is  a  function  of  temperature,  humidity,  and  air  pressure, 
where  within  the  PBL,  the  effects  of  air  pressure  can  be  ignored.  Refractive  index  fluctuations 
can  be  expressed  as  a  linear  combination  of  fluctuation  in  temperature  and  humidity.  Equations 
(3.1)  and  (3.2)  were  used  to  obtain  the  scalar  structure  function  parameters  expressed  as 

C2  =  (6xl04Pn2)2C2a2  (3.6) 

where  Cq2  is  the  humidity  structure  function  parameter  (Wesely,  1976)  and 

ar2  =  l-  2arTqC/Cq  +  a2Cr2/C2  and  (3.7) 
rr,  =  <V(CyC?).  (3.8) 

In  these  expressions,  the  quantity  (a2  - 1)  gives  the  deviations  from  the  pure  moisture  dominated 
6>j2.  Cj  is  the  temperature  structure  function  parameter ,  rTq  is  the  temperature-moisture 
correlation  coefficient,  and 

a  =  0.13(1  +  15.4q/T)  (3.9) 

(Fairall  1991).  Contributions  of  CT2  and  C2  to  C2  were  investigated  by  Wesely  (1976).  Within 
the  lower  troposphere,  he  found  C2  dominated.  VanZandt  (1978)  also  extended  this  result 
above  the  boundary  layer  to  a  few  kilometers  of  the  troposphere  based  on  a  rawinsonde,  radar, 
and  model  data  correlation  study.  Above  a  few  kilometers  and  into  the  ionosphere,  C2  is  the 
primary  contributor  to  C2.  Exceptions  within  the  PBL  for  C2  domination  may  exist  in 
extremely  dry  and  strongly  convective  atmospheres.  Therefore,  a2  =  1  for  most  marine  and 
continental  PBLs  in  the  tropics  and  mid-latitudes. 

Following  White  et.  al  (1991),  after  rearranging  (2.6)  in  terms  of  C2,  substituting  (2.5) 
into  (2.6),  and  where  T=  268°,  P*P0ef'm  is  the  pressure  in  millibars,  z  is  the  altitude  in 
kilometers,  P0  =  1000  mb,  and  H  =  8  km  yields 

C2  =  (1.66 x  103r/P)2a;2(Tj/0.38)AI/3.  (3.10) 
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From  this  equation,  the  turbulent  humidity  structure  of  the  atmosphere  can  be  calculated. 

In  order  to  relate  t]  to  a  quantity  measured  by  the  Doppler  radar  we  use  the  reflectivity  equation 
given  by  VanZandt  et  al.  (1978): 

n  =  9nck0B(aTc  +  TJ(R2/AR2)(SNR)[2a2PlFrAe  cos(x)]'1  (3.11) 
where  c  =  2.998  x  10s  m  s"1  is  the  speed  of  light  in  a  vacuum,  k0  =  1.3803  x  10"23  J  mol"1  K'1  is 
Boltzmann's  constant,  R  is  the  range,  and  the  remaining  parameters  are  listed  in  Table  3. 


cm 

Parameter 

Value 

K 

Effective  antenna  area  (one  section  equal  to  geometric  area 

47.0  m2 

a 

Efficiency  of  radar  antenna  and  transmission  line 

0.10 

p, 

Peak  pulse  power 

9.0  kw 

X 

Off  vertical  beam  axis  angle 

14.5° 

T 

*  rx 

Receiver  noise  temperature 

100  K 

m 

Cosmic  noise  temperature 

10-100  K 

AR 

Range  resolution 

149.9  m 

Pr 

Pulse  repetition  frequency 

10kHz 

B 

Bandwidth  of  integrating  filter 

50.0  Hz 

Table  3.  Radar  characteristics  and  signal  processing  parameters  used  in  deriving  the  relationship  between 
reflectivity  and  signal  to  noise  ratio  (SNR)  for  the  Penn  State  404  MHz  wind  profiler.  (Penc,  1995) 

Combining  (2. 10)  with  (2. 1 1)  and  substituting  values  from  Table  3  gives  the  form  of  the  C2 

equation  used  in  this  study: 

C2  =  1.63  x  1 0',s(T2/P)2  ar'2R2l  [fNmo  (2.12) 

where  SNR  is  now  conveniently  in  decibels  (dB),  as  it  was  recorded  during  the  four  LOWS  IOPs 
used  in  this  study. 

4.  Data  Analysis  Procedures 

A  number  of  steps  were  taken  organize  the  raw  radar  and  rawinsonde  data  into  a  format 
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suitable  for  computational  and  visual  analysis.  In  order  to  compute  entrainment,  C2  values  had 
to  be  derived  from  the  raw  radar  data.  To  compute  moisture/temperature  contributions  to  C2, 
a2  must  be  derived  from  rawinsonde  data.  For  convenience  sake,  personal  computer-based 
display  and  spreadsheet  software  were  used. 

4.1  Radar  Data 

Originally,  the  raw  radar  data  from  the  North  Rose  wind  profiler  was  stored  on 
magnetic  media  as  individual  files  sorted  by  day.  Most  of  the  files  contained  an  entire  day's  data 
(>  4  Mb)  with  information  from  each  beam  (east,  north,  and  vertical)  and  mode  (low,  high)  that 
was  recorded  every  three  minutes  (sub-hourly  data).  The  sub-hourly  data  incorporated  raw 
spectra  and  moments  tables,  consisting  of  the  following  parameters  in  tabular  form  by  range 
gate:  horizontal  velocity,  returned  power,  noise,  signal-to-noise  ratio  (SNR),  and  spectral  width. 
In  this  boundary  layer  growth  study,  only  the  SNR  fields  from  the  moments  tables  for  the  vertical 
beam  were  required  for  computations  and  analysis.  Additionally,  the  higher  resolution  low  mode 
data  was  desired  as  the  boundary  layers  in  the  study  remained  below  the  low  mode  top  range  gate 
(2600m).  Since  only  a  small  portion  of  the  copious  raw  radar  data  was  required,  a  C  program 
was  written  to  parse  and  store  the  applicable  data  on  a  desktop  PC.  Also,  the  output  of  the  DOS- 
based  executable  needed  to  be  in  a  grid  format  compatible  with  SURFER,  a  mathematical  grid 
manipulation  and  isoplething  software  package. 

Equation  (3.12)  was  applied  to  the  original  SNR  ratio  grids  to  obtain  C2  grids,  from 
which  Iog(C?2)  grids  were  calculated.  Next  a  2, 1  time-height  matrix  smoothing  was  applied  to 
the  log(C?2)  grids  and  then  plotted. 

4.2  Rawinsonde  Data 

Raw  rawinsonde  data  was  used  to  determine  temperature  and  moisture  contributions  to 
C2  and  combined  with  log(C?2 )  values  to  obtain  entrainment  velocities  (we).  The  original  ASCII 
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files  were  imported  into  MS  EXCEL  spreadsheets,  then  meteorological  parameters  were  derived. 
To  make  time-height  series  spreadsheets,  like  parameters  from  individual  soundings  were 
combined.  Sparse  linear  data  interpolation  was  necessary  to  maintain  consistent  height  data, 
however  because  the  resolution  of  the  sounding  data  (~50m)  is  significantly  better  than  the  radar 
data  in  low  mode  (~150m),  there  should  be  little  error  injected  into  the  entrainment 
computations. 

The  meteorological  parameters  were  calculated  from  the  original  rawinsonde  data  as 
shown  in  T able  4,  where  RH  is  relative  humidity,  p  is  pressure  (mb),  pa  is  surface  pressure  (mb), 
T  is  temperature  ( °Q\ 


specific  humidity  ( q ) 

q  =.  037990516  RHexp[  (17.27  T)/(T  +  273)  -  35.86)] /p 

potential 
temperature  (0) 

0=(T  +  273)(p/p)m 

saturation  vapor 
pressure  (es) 

es  =  6.1078  exp [  (17.27  T)/(T  + 273)  -  35.86)] 

saturation  mixing 
ratio  (wt) 

w,  =  0.622  (e/(p  -  ej) 

actual  mixing 
ratio  (w) 

w  =  w,  RH/ 100 

virtual  potential 
temperature  (6J 

ev  =  d(i  +  o.6i  w) 

Table  4.  The  meteorological  parameters  calculated  from  the  original  rawinsonde  data  and  applicable 
equations. 

4.3  Alpha  ( ar2)  Calculations 


To  accurately  assess  the  individual  temperature  and  humidity  structure  parameter 
contributions,  a2  must  be  estimated  using  (3.12),  as  the  quantity  (a2  -1)  gives  the  deviation 
from  pure  moisture  dominated  C2.  It  is  generally  agreed  upon  that  a2  is  near  unity  in  the  non- 
convective  PBL.  However  in  the  LOWS  cases  studied,  the  PBL  is  highly  convective,  therefore  a 
detailed  analysis  was  required  using  sounding  data.  From  Burk  (1981),  we  have  a  method  for 
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determining  relative  contribution  to  C2  using  the  following  equation,  which  is  another  version  of 
(3.7): 

ar2  =  1  -  [(2Ae/C)  +  (4ee/pT)][2.2  -  2.4Tr/3.9pr]  +  [( Ae/C )  +  (4. 1) 

(2ee/pT)]2[(0.5  -  2.67V  +  1.4(Tr)2)/(3.9(pr)2)] 

where  A  is  a  coefficient  in  microwave  refractivity  (77.6  x  10'6  K  mb'15,  C  is  a  coefficient  in 

microwave  refractivity  (0.375  K2  mb'1),  p  is  pressure  (mb),  T  is  absolute  temperature,  e  is  vapor 

pressure  (mb),  eis  the  molecular  weight  of  water/dry  air  (0.622),  p  is  the  density  of  air,  r  = 

Aq/ A  0V  across  the  interfacial  layer.  To  obtain  the  changes  of  specific  humidity  and  virtual 

potential  temperature  at  a  specific  height,  gradients  were  calculated  from  rawinsonde  data  using 

respective  values  immediately  below  and  above  the  layer  in  question.  All  computations  were 

accomplished  on  spreadsheets  from  which  time-height  series  graphs  were  created,  then  imported 

to  SURFER  for  analysis. 

4.4  Entrainment  Velocity  (wO  Calculations 

Entrainment  is  the  mixing  down  of  less  turbulent  air  into  the  PBL,  thereby  enhancing 

boundary  layer  growth.  (Stull,  1988)  The  entrainment  velocity  is  therefore  one  of  the  key 

parameters  to  evaluate  interactions  occurring  in  the  inversion  near  the  cloud  tops.  After  applying 

the  moisture  deviation  from  (4.1)  to  C2  values  obtained  from  (3.12),  entrainment  values  can  be 

calculated  following  Wyngaard  and  LeMone  (1980): 

we  =  6  A0V  s1'3  C2  [11.2  ( Aq)2!]'1  (4.2) 

evaluated  at  the  middle  of  the  interfacial  layer  where  0V  is  the  virtual  temperature,  /'is  the  lapse 
rate  of  0y  above  the  inversion,  A  denotes  changes  across  inversion,  e  is  the  rate  of  dissipation  of 
turbulent  kinematic  energy  for  cloudy  conditions  (1.0  x  10'3m2s'3).  The  base  and  the  top  of  the 
inversion  ware  taken  to  be  the  span  of  the  lowest  temperature  inversion  identified  in  virtual 
temperature  profiles  and/or  visual  analysis  of  moisture  profiles.  /  was  calculated  for  the  first 
1000  m  above  the  inversion. 
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5.  Interpretation  of  Data 


Entrainment  of  air  from  the  free  atmosphere  into  the  convective  boundary  layer  enables 
larger  precipitation-producing  convective  eddies  in  the  mixed  layer.  In  order  to  obtain  reliable 
estimations  of  the  height  of  the  convective  boundary  layer  and  entrainment  velocitiy  from  the 
radar  data,  accurate  moisture  profiles  were  needed.  According  to  the  theory  presented  in  section 
2,  the  radar  refractivity  is  proportional  to  gradients  of  moisture  and  temperature.  It  has 
previously  been  determined  that  in  Type  II  CTBLs,  the  radar  reflectivity  of  clear  air  is  primarily 
dependent  upon  moisture  contributions  (Fairall,  1991).  The  humidity  structure  function 
parameter  serves  as  a  relative  measure  of  these  gradients  of  turbulent  mixed  moisture  and  can  be 
calculated  from  radar  reflectivity  data.  Since  it  is  known  that  moisture  generally  decreases  with 
height  and  the  moisture  gradients  are  largest  at  the  top  of  the  CTBL,  time-height  graphs  of 
log(C?2)  values  can  be  used  to  monitor  CTBL  depths.  In  order  to  ensure  this  relationship  is 
valid,  either  a  parameterization  concerning  the  contributions  of  moisture  to  the  index  of 
refraction  must  be  made  or  a  calculation  from  rawinsonde  data.  The  radar  humidity  data  has  the 
temporal  advantage  of  observations  every  3  1/2  minutes,  whereas  rawinsondes  were  launched  at 
various  frequencies  (1  to  5  per  day).  However,  in  the  vertical,  the  rawinsonde's  humidity  data 
resolution  was  50  m,  compared  with  150  m  for  the  radar. 

Signal-to-noise  ratio  (SNR)  time-height  series  data  appears  in  Figs.  1-4  for  the  11  -14 
January  single  band  events  and  in  Figs.  5-8  for  the  25-26  February  and  28  February  - 1  March 
muliple  band  events.  Figs.  1-8  show  values  of  SNR  generally  decrease  with  height,  with  a  few 
exceptions.  One  major  exception  is  in  Fig  3  for  13  January,  where  the  radar  malfunctioned  for  a 
short  period  around  1900  UTC.  Without  the  known  temperature  and  moisture  contributions  to 
the  refractive  index,  other  direct  interpretations  as  related  to  convective  boundary  layer 
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features  are  premature.  Therefore,  the  in-depth  discussion  of  radar  data  will  be  later  in  this 
section. 

5.1  Alpha  <  a.2)  Calculations 

As  previously  stated,  the  quantity  (a2  - 1)  gives  the  deviation  from  pure  moisture 
dominated  C2.  Studies  of  non-convective  PBLs  have  shown  a2  near  unity.  Model  studies  have 
produced  values  for  a2  near  unity  for  cloud-free  convective  boundary  layers.  According  to 
Fairall  (1990)  there  are  three  ways  clouds  change  the  radar  refractivity  potential  of  the  boundary 
layer:  an  additional  thermodynamic  coupling  between  temperature  and  moisture  results  due  to 
condensation  processes,  the  entrainment  rate  is  greatly  enhanced  by  cloud-top  radiative  cooling, 
and  cloud  liquid  water  droplets  provide  an  additional  scattering  mechanism  for  radar.  Fairall 
states,  based  on  simple  physical  reasoning,  the  results  for  a2  are  approximately  valid  for  the 
entrainment-induced  effects  associated  with  clouds.  Similar  conclusions  were  reached  in  a 
modeling  study  of  two  maritime  boundary  layer  cases  even  in  the  presence  of  clouds  by  Burk 
(1981). 

The  results  of  a2  for  the  LOWS  data  are  best  shown  by  Fig.  9.  Throughout  the  time 
frame,  values  for  a2  ranged  from  1.02  to  1.04  within  the  boundary  layer.  Therefore,  based  on 
(3.6),  contributions  of  C2  to  C2 are  96  to  98  percent.  Additionally,  for  modeling  studies  and 
calculations  involving  the  refractive  index  structure  parameter,  an  assumption  of  a2  near  unity 
will  introduce  very  little  error  (<  4%).  Of  note  in  Fig.  9,  as  the  boundary  layer  depth  decreased 
under  the  effects  of  subsidence,  a2  showed  increasing  variability  above  the  boundary  layer.  By 
visual  inspection,  this  result  appears  to  be  directly  correlated  with  the  specific  humidity  vertical 
profile  in  Fig  10.  However,  these  moisture  gradients  still  have  little  effect  on  the  overall 
contributions  of  C2  to  C2  within  the  convective  boundary  layer. 
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5.2  Humidity  Data 


A  visual  comparison  of  specific  humidity  (Fig.  10)  with  log(C42)  for  1 1  - 14  January 
(Figs.  11-14)  shows  excellent  correlation.  Specifically,  Fig.  10  at  2300  UTC  on  1 1  January 
shows  increased  specific  humidity  at  600m,  1600m,  and  2500m.  Whereas,  Fig.  1 1  shows 
increased  log (C/)  values  at  600m,  1600m,  and  2200m.  This  corresponded  with  the  passage  of 
a  surface  front,  after  which  lake-effect  snow  bands  formed  in  the  post-frontal  environment.  Cold 
front/trough  passages  were  typified  by  moderately  narrow  uniform  vertical  log(C42))  fields.  In 
Fig.  11, 1900-2100  UTC  on  11  January,  uniform  vertical  values  of  log(C42)  exceed  -3.6  then 
rapidly  transitioned  to  -4.5.  Light  snowfall  in  the  pre-frontal  environment  was  seen  as  high 
values  (>  -4.2)  of  log(C?2)  with  loose  vertical  gradients  and  a  boundary  layer  depth  of  around 
2200m.  Three  distinct  snow  bands  are  seen  in  Fig.  12  on  12  January,  0400-0800  UTC  indicated 
by  very  narrow  uniform  vertical  log(C42))  fields. 

Continuing  the  visual  comparison,  Fig.  10  at  1600  UTC  on  13  January  indicates 
specific  humidity  increases  at  1000-1200m,  1600-2000m,  and  2250-2450m.  Whereas,  Fig.  13 
shows  increased  log(C42)  values  1000-1300m,  1600-1 800m,  and  2200-2400m.  In  Fig.  14  on  14 
January,  the  last  snow  band  is  seen  around  0900  UTC,  then  frontal  passage  at  1400  UTC,  with 
considerable  drying  of  the  atmospher  thereafter.  Fig.  15  depicts  specific  humidity  for  the  28 
February,  1700  UTC  sounding.  While  the  change  in  specific  humidity  is  nearly  linear,  closer 
inspection  reveals  decreases  in  the  rate  of  change  at  500-800m,  1250-1600m,  2250-2450m.  Fig 
16  shows  a  similar  pattern,  however  visual  correlation  is  not  as  prominent  as  in  the  previous 
examples. 

5.3  Boundary  Laver  Heights 

The  depth  of  the  boundary  layer  limits  convective  eddys  and  therefore  the  intensity  of 
lake-effect  snow  events.  From  vertically  pointed  radars  we  can  obtain  moisture  profile  (time- 
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height  series  log(C?2))  data  that  may  be  used  to  estimate  boundary  layer  heights.  Normally, 
atmospheric  moisture  content,  and  hence,  C2  decreases  with  height.  However,  sharp  moisture 
gradients  near  the  inversion  cause  local  C 2  maximiums,  from  which  depth  of  the  boundary  layer 
can  be  distinguished.  For  the  most  part,  these  features  are  easily  recognized.  As  one  might 
expect,  the  banded  snow  events  that  frequent  the  North  Rose  site  can  be  seen  as  sharp  changes  in 
boundary  layer  depth. 

The  25-26  February  multiple  snow  band  event  exemplifies  the  approach  to  determining 
boundary  layer  heights.  In  Fig.  19  at  around  0200  UTC,  a  rapidly  moving  "Alberta  Clipper"  was 
transitioning  through  the  area  producing  widespread  light  snow.  The  boundary  layer  remained 
around  2300m  until  approximately  0800  UTC  when  cold  dry  air  moved  rapidly  in  behind  the 
system  as  it  exited.  By  1200  UTC,  the  boundary  layer  depth  was  reduced  to  1200m.  Despite 
the  short  fetch  due  to  northerly  surface  winds,  convection  was  seen  around  1630  UTC  as  the 
boundary  layer  quickly  rose  to  2300m  when  a  lake-effect  snow  band  moved  over  the  radar  site. 
At  2100  UTC  the  boundary  layer  dropped  back  to  1200m.  Then  the  last  snow  band  passed 
overhead  0230-0400  UTC,  after  which  the  boundary  layer  rapidly  deflated  as  strong  high 
pressure  and  subsidence  moved  in.  This  event  left  4  cm  of  snow  at  the  radar  site. 

5.4  Entrainment  Velocity  (V.) 

The  entrainment  velocities  in  Table  5  are  within  a  factor  of  2  -  3  of  calculations  made  by  Penc 
(1995)  for  the  LOWS  data.  In  the  Penc  entrainment  velocity  calculations,  a  lake-scale  average 
for  each  major  snow  event  was  obtained  using  a  kinematic  equation  for  growth  of  a  convective 
boundary  layer  initially  presented  by  Deardorff  and  Peterson  (1980).  This  technique  uses 
boundary  layer  height  differences  averaged  across  the  lake  to  compute  growth  rates,  from  which 
entrainment  is  derived.  Our  calculations  apply  only  to  the  North  Rose  radar  site,  but  are  within 
expected  ranges. 
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Results  of  Entrainment  Velocity  (wr1  Calculations 
12-13  Jan  90 
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Table  5.  Values  of  variables  derived  from  radar  and  rawinsonde  data  from  12-13  January  1990  used  in 
calculating  entrainment  velocities  along  with  their  units. 

The  entrainment  velocities  in  Table  5  are  within  a  factor  of  2  -  3  of  calculations  made  by  Penc 
(1995)  for  the  LOWS  data.  In  the  Penc  entrainment  velocity  calculations,  a  lake-scale  average 
for  each  major  snow  event  was  obtained  using  a  kinematic  equation  for  growth  of  a  convective 
boundary  layer  initially  presented  by  Deardorff  and  Peterson  (1980).  This  technique  uses 
boundary  layer  height  differences  averaged  across  the  lake  to  compute  growth  rates,  from  which 
entrainment  is  derived.  Our  calculations  apply  only  to  the  North  Rose  radar  site,  but  are  within 
expected  ranges. 

Table  6  shows  calculated  entrainment  velocities  for  two  multiple  snow  band  events. 
Due  to  limited  rawinsonde  data  for  the  25-26  February  and  28  February  - 1  March  1990  multiple 
snow  band  events,  entrainment  velocities  could  only  be  calculated  for  two  time  periods.  The  25 
February  at  0300  UTC  we  calculation  result  appears  rather  large 
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Results  of  Entrainment  Velocity  (wc)  Calculations 
25  Feb  90  and  28  Feb  90 
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Table  6.  Values  of  variables  derived  from  radar  and  rawinsonde  data  from  25  February  and  28  February 
1990  used  in  calculating  entrainment  velocities  along  with  their  units. 

compared  with  the  results  from  Table  5.  However,  a  rough  order  of  magnitude  calculation 
derived  from  Fig.  19  between  0200  and  0300  UTC  reveals  the  boundary  layer  grew 
approximately  400m,  or  the  entrainment  velocity  averaged  1 1.2  cm/sec.  Because  of  the 
consistency  between  wt  determined  from  two  approaches,  the  instantaneous  entrainment  velocity 
result  for  0300  UTC  appears  to  be  reasonable. 

6.  Conclusions 

Consistent  with  previous  findings  by  White  et  al.  (1991)  for  Type  II  CTBLs  and 
cloudless  convective  boundary  layers,  gradients  of  moisture  are  the  primary  contributor  to  the 
radar  index  of  refraction  within  a  Type  I CTBL,  and  therefore  can  be  reliably  sensed  using 
vertical  beam  radars.  The  increased  temporal  resolution  of  the  radar  moisture  analysis  makes 
this  an  excellent  tool  for  monitoring  cloud-topped  boundary  layer  depths.  Moisture  profiles 
obtained  from  satellites  are  sufficient  for  sensing  the  upper  and  mid-levels  of  the  atmosphere. 
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However,  because  of  the  poor  resolution  of  satellite-based  moisture  sensors  for  boundary  layer 
application,  wind  profiler  networks  will  remain  necessary.  Boundary  layer  height  and  cold 
front/trough  passage  can  be  gauged  using  the  humidity  data  from  these  wind  profiler  radars. 
Additionally,  using  the  assumption  (a2  -  1)  for  a  Type  I CTBL  is  valid  and  will  introduce  very 
little  error  ( <  4%)  in  radar  applications  and  modeling  studies.  Entrainment  velocity  rates  and 
boundary  layer  growth  can  be  reliably  estimated  using  humidity-derived  SNR  radar  data,  as 
entrainment  velocity  calculations  for  the  North  Rose  site  were  consistent  with  a  lake-averaged 
kinematic  equation  for  growth  of  the  convective  boundary  layer.  Also,  rough  order  of  magnitude 
entrainment  velocity  rates  obtained  from  time-averaged  radar  data  appeared  to  validate  the 
instantaneous  we  calculations. 
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